Fine-scale oscillations in the wavelength- and intensity-dependence of high-harmonic 

generation: connection with channel closings 
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We investigate the connection of recently identified fine-scale oscillations in the dependence of the 
yield of the high-harmonic generation (HHG) on wavelength A of a few-cycle laser pulse [K. Schiessl, 
K.L. Ishikawa, E. Persson, and J. Burgdorfer, Phys. Rev. Lett. 99, 253903 (2007)] to the well- 
known channel closing (CC) effect. Using the Lewenstein model of HHG, we identify the origin of the 
oscillations as quantum interference of many rescattering trajectories. By studying the simultaneous 
variations with intensity and wavelength, different models for the interference of channel closing 
peaks can be tested. Contrary to theoretical predictions for short-range potentials, the peaks are 
not located neither at nor just below the CC condition, but a significant shift is observed. The long 
Coulomb tail of the atomic potential is identified as the origin of the shift. 

PACS numbers: 32.80.Rm,42.65.Ky,32.80.Fb 
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I. INTRODUCTION 



High harmonic generation (HHG) represents a versatile 
and highly successful avenue towards an ultrashort co- 
herent light source covering a wavelength range from the 
vacuum ultraviolet to the soft X-ray region [l[. This 
development has opened new research areas such as at- 
tosecond science [2 H[ and nonlinear optics in the XUV 
region The fundamental wavelength A used in most 

HHG experiments to date is in the near-visible range 
(~ 800 nm). The cutoff law for the harmonic spectrum 
Ip + 3.17t/ p , where I p denotes the binding energy 

2 /A , ,2 _ ^2x2/^2 the 



E c 

of the target atom and U p = Fq/4u> 2 = F Q I \ 2 /16ir 
ponderomotive energy (Fq: laser electric field strength), 
suggests that a longer fundamental wavelength would be 
advantageous to extend the cutoff to higher photon en- 
ergies, since U p increases quadratically with A. This has 
stimulated an increasing interest in the development of 
high-power mid- infrared (~ 2 (im) laser systems, e.g., 
based on optical parametric chirped pulse amplification. 
The first generation of water-window harmonics with 
clear plateau and cut-off structures has recently been 
reported [9(. Along those lines the dependence of the 
HHG yield on A has become an issue of major interest 
S 0, H E3, [3 • It had long been believed that the spread- 
ing of the returning wavepacket would result in a A -3 
dependence of the HHG efficienc y fill as long as ground 
state dep letion can be neglected [I4I; experimental find- 
ings [13[ provided partial support. Recently, however, 
Tate et al. [T3| have reported a different wavelength- 
scaling of HHG between 800 nm and 2 /im calculated 



with the time-dependent Schrodinger equation (TDSE) 
for Ar and a strong-field approximation (SFA) for He. 
They found the yield to be described by a power-law 
oc \~ x with 5 < x < 6. Investigating the A depen- 
dence on the level of single-atom response for H and Ar 
by numerically solving the time-dependent Schrodinger 
equation we could confirm the overall scaling with an in- 
verse power law exceeding five @ , the harmonic yield was 
found not to depend smoothly on the fundamental wave- 
length, but to exhibit surprisingly rapid oscillations with 
a period of 6 — 20 nm depending on the wavelength region. 
A semiclassical analysis based on the SFA has revealed 
that the rapid oscillations are due to the interference of 
five to ten different rescattering trajectories Q. More- 
over, we found the oscillations to be stable with respect 
to variations of the pulse envelope as long as the effective 
pulse length and thus the number of relevant trajectories 
remains equal, while the amplitude of the oscillations de- 
creases with decreasing pulse length These observa- 
tions underscored the view that the oscillations are due 
to the interference of quantum paths. 

Oscillations of the HHG yield have previously been re- 
ported in terms of the dependence on the intensity of the 
driving laser in oc Fq , both experimentally fl5l . Il6| and 
theoretically [13, EH • Borca et al. and Milosevic and 
Becker [l8| have shown that HHG is enhanced at channel 
closings (CC), i.e., when 



R= I * 



(i) 
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is an integer. Channel closing in this context refers to 
the threshold for multiphoton ionization in a laser field. 
Most of these theoretical studies employed zero-range po- 
tentials or the SFA which both neglect the influence of 
the long-range potential on the ionized electron. 
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Frolov et al. Q have recently analyzed the wavelength- 
dependence of HHG in terms of channel closings (or 
threshold phenomena). They have calculated the har- 
monic yield using the time-dependent effective range the- 
ory, and shown that the peaks of the yield oscillation 
around A = 1 /im coincide with integer values of R if 
an effective ionization potential I p (e.g., 10.5 eV for H) 
is used in place of I p in Eq. fTJ). This method is, how- 
ever, strictly applicable only for short-range potentials 
and also neglects the excited atomic states. On the other 
hand, we have recently found [1| channel closing peaks 
in the TDSE-calculated HHG yields around 1 /j,m and 2 
/Ltm wavelengths which are characterized by a spacing of 
5R = 1, as expected from the CC picture, when the true 
ionization potential I p is used. 

In the present paper, we study the connection between 
the oscillation in the wavelength-dependence of the HHG 
yield and the channel-closing in more detail. We study 
the harmonic spectrum in the two-dimensional parame- 
ter space of intensity I and driver wavelength A c . We 
compare the results of the full 3D TDSE solution with 
the strong-field approximation and a truncated Coulomb- 
potential model in order to delineate underlying mecha- 
nisms. We find that the correspondence of the modula- 
tion period to SR = 1 holds for a wide wavelength range 
between 800 nm and 2/im, and that the peak positions in 
terms of R are almost independent of laser intensity. The 
systematic displacement of the peak positions relative to 
integer values is found to be consistent with the effect of 
the long-range Coulomb tail on the returning electron. 

The present paper is organized as follows. Section |TT] 
summarizes the two complementary integration schemes 
employed for a full numerical solution of the TDSE. In 
Sec. IIIII we discuss the overall wavelength dependence at 
a fixed value of fundamental intensity all the way from 
A =800 nm to 2/im. We also analyze small-scale oscilla- 
tions in terms of quantum-path interference based on the 
saddle-point analysis (SPA) [H QJ]. In Sec. |TV]we in- 
vestigate the variation of the A-dependence of HHG with 
intensity and pulse shape. In Sec.|V]we discuss the period 
of the oscillations in terms of the channel-closing num- 
ber R and investigate its robustness against the variation 
of the wavelength region, the driver intensity, and pulse 
shape. In Sec. I VII we discuss the origin of the peak shift 
from integer R values and clarify how the Coulomb tail 
of the atomic potential affects the peak position. Con- 
clusions are given in Sec. IVIII Atomic units are used 
throughout the paper unless otherwise stated. 



II. NUMERICAL METHODS 

We solve the atomic time-dependent Schrodinger equa- 
tion (TDSE) in the length gauge for a linearly polarized 
laser field with the central wavelength A c = 2ttc/uj, 



where F(t) = Fo f (t) sm(u)t) denotes the laser electric 
field, f(t) is the envelope function and V e g(r) the atomic 
potential. For hydrogen (H), V e g(r) is the bare Coulomb 
potential while for argon (Ar) we employ a model poten- 
tial [19( within the single-active electron approximation 
which reproduces the binding energy to an accuracy of 
typically w 10~ 3 . We employ two complementary meth- 
ods to solve Eq. @ in order to establish reliable and 
consistent results. 

In the first method, Eq. @ is numerically integrated 
using the alterna ting direction implicit (Peaceman- 
Rachford) method [20| with a uniform grid spacing Ar 
being dependent on the numerical problem in the range 
of 10~ 2 < Ar < 6.25 x 10~ 2 a.u. In general, a finer grid 
spacing is needed for a longer wavelength, and also for 
Ar than for H. In order to reduce the difference between 
the discretized and analytical wave function, we scale the 
Coulomb potential by a few percent at the first grid point 
[Ij. The time step At is typically 1/16000 of an opti- 
cal cycle for 800 nm wavelength, i.e., 6.895 x 10~ 3 a.u.. 
This algorithm is accurate to the order of 0(At 3 ). In 
the second method, the TDSE is integrated on a finite 
grid by means of the pseudo-spectral method [22| which 
is also accurate to the order of 0(Ai 3 ). It allows for 
larger time steps of the order of 0.1 atomic units. The r 
coordinate is discretized within the interval [0, r max ] with 
a non-uniform mesh point distribution. The innermost 
grid point is typically as small as 2.5 x 10~ 4 a.u., enabling 
an accurate description near the nucleus. A smooth cut- 
off function is multiplied at each time-step to avoid spu- 



rious reflections at the border 



while another cut- 



off function prevents reflections at the largest resolved 
energy -E max - For Ar the occupied states supported by 
the model potential are dynamically blocked during the 
time evolution by assigning a phase corresponding to an 
unphysically large and positive energy eigenvalue [23|. 
We calculate the dipole acceleration d(t) — — <9 2 (z(i)), 
employing the Ehrenfest theorem throu gh the relation 
d\t) = (ip(r,t) | cos6>/r 2 - F(t) | ^(r,t))Q, where the 
second term can be dropped as it does not contribute to 
the HHG spectrum. 

For the wavelength-dependence of the harmonic yield, 
in particular the global scaling, it is important to specify 
the definition of the integral yield. One can focus on a 
given number of harmonic orders, on a fixed energy inter- 
val, or the entire spectrum. Following Refs. @, 0, H, [3 , 
we consider in this work the HHG yield defined as energy 
radiated from the target atom (single-atom response) per 
unit time [25| integrated for a fixed photon energy range, 
specifically from 20 to 50 eV: 



Ay = 



1 



3c 3 T 



50 cV 



|a(w)| 2 du> , 



(3) 



20 cV 



.^(r,t) = 



|v 2 +V eS (r)+zF(t) 



V-M), (2) 



where T denotes the pulse duration. Note that the energy 
window AE of the output radiation (here 20 to 50 eV) 
is kept constant when analyzing Ay as a function of 
A. Clearly, both the number and order of the harmonic 
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FIG. 1: (Color online) Wavelength dependence of the inte- 
grated harmonic yield AY between 20 and 50 eV as a func- 
tion of (a) wavelength A and (b) channel-closing number R. 
■ : TDSE results obtained by the pseudo-spectral method. 
In addition, results on a fine scale (solid line) are presented. 
Dashed line: fit AY oc A -5 ' 3 . 



peaks lying in the fixed energy interval change as A c is 
varied. 



III. WAVELENGTH DEPENDENCE 

We adopt the laser parameters of Ref . [l4[ , with a fixed 
peak intensity of 1.6 x 10 14 W/cm 2 , a variation of A be- 
tween 800 nm and 2 /xm, and an envelope function f(t) 
corresponding to an 8-cycle flat-top sine pulse with a half- 
cycle turn-on and turn-off. 

Figure Q] displays the HHG yield for atomic hydrogen 
calculated on a fine mesh in A with a spacing of 1 nm. 
Superimposed on a global power-law dependence AY oc 
\~ x (i « 5) | 0,11, EJ], we find remarkably strong 
and rapid fluctuations through the entire A range. The 
origin of this oscillation can be identified as the quantum 
interference of up to ten rescattering trajectories, based 
on the SFA analysis @,|||. 

The effect of the interference of multiple quantum 
paths was previously studied in the context of the inten- 



sity dependence of HHG and above-threshold ionization 
(ATI) [16|, EE El| • Using the quasistationary quasienergy 
state theory for a zero-range potential and the strong- 
field approximation, Borca et al. 117], and Milosevic and 
Becker [18[ have shown that HHG exhibits resonance-like 
enhancement when A-photon ionization channel is closed 
with increasing intensity, i.e. the parameter R (Eq. (fTj)) 
becomes an integer. Zair et al. [161 ] have very recently 
reported experimental observation of the oscillation in 
the intensity dependence of the HHG yield as evidence 
of the interference between the short and long paths. The 
analysis in Ref. may be applied to the wavelength- 
dependence as well, suggesting to look at our results in 
terms of R. Data of Fig. [lja) are replotted in terms of 
R in Fig. HJb), which permits a detailed analysis of the 
channel closing behavior (see below). 

In Figs. [2] and [3] we reexamine the role of quantum 
paths in the oscillations of the wavelength-dependence of 
the harmonic yield as a function of R on a finer R scale. 
We compare full TDSE solutions with approximations 
based on the strong field approximation (SFA) [ll|, [HI- 
We first apply the Gaussian model |Tlj], in which the 
ground-state wave function has the form 



*<"-(;) 



3/4 



-ar 2 /2 



(4) 



where a is chosen to reproduce I p . An appealing point of 
the Gaussian model is that the dipole transition matrix 
element also takes a Gaussian form 111, 



d(p) = i 



1 



3/4 



£ e -p 2 /2c 



(5) 



and that one can evaluate the integral with respect to mo- 
mentum in the formula for the dipole moment (Eq. (8) 
of Ref. [ll[) analytically, without explicitly invoking the 
notion of quantum paths. Unphysically rapid decrease 
for p 2 /2a ^> 1 limits the application of Eq. ([5]) to har- 
monic orders with momenta of the returning electron not 
substantially exceeding p 2 /2cv « 1. We have confirmed 
that the resulting harmonic spectra have an adequate 
plateau and cut-off structure for the value of a (1 and 
2 a.u.) used in the present study. The obtained wave- 
length dependence of the HHG yield, expressed in terms 
of R (Fig. [2] (b)), exhibits oscillations similar to that in 
the TDSE result (Fig. [2] (a)), although peaks are found - 
contrary to TDSE results - near integer values of R. In 
addition we employ complex solutions of the saddle-point 
approximation (SPA) [18| . while we have previously ob- 
tained similar results by employing classical trajectories 
0, [1|. Up to 16 possible trajectories for each individ- 
ual photon energy are considered. When including up 
to ten and twelve returning paths for the case of H and 
Ar, respectively, the SPA can reproduce the modulation 
depth and frequency of the A oscillations of the TDSE 
and the Gaussian model reasonably well, thus strongly 
supporting the quantum path interference as the origin 
of the fluctuations. The SPA result for Ar with twelve 
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trajectories (Fig. [3] (b)) reproduces even the small peaks 
between the main peaks. 

We emphasize the remarkable variation on a fine A 
scale. One might suspect that the oscillation as in Figs. 
[JH3] would be specific to monochromatic driver pulses 
and smeared out for the case of ultrashort broadband 
pulses. The pulse shape used in this study is, however, 
not monochromatic but its spectral width AA is ~ 10% 
of A c . The rapid variations of the harmonic yield occur 
on a scale SX much smaller than this width. This find- 
ing, at a first glance surprising, is a direct consequence 
of the quantum path interference. It follows from the ex- 
istence and the fixed spacing in between discrete points 
in time - controlled by A c - at which electronic trajecto- 
ries are launched. As long as the few-cycle pulse permits 
the generation of a set of a few quantum paths in subse- 
quent half-cycles, the overall temporal characteristics of 
the driver pulse is of minor importance. We have also 
checked that the fluctuations in the harmonic yield are 
not an artefact of our particular choice of f(t). They 
can be observed also for "smoother" pulse shape such 
as sin 2 and Gaussian pulses as well as shorter pulses, 
provided that the pulse can support multiple returning 
trajectories (Fig. [3]). The temporal profile of the pulse 
influences the detailed shape of the interference pattern, 
in particular the amplitude of the oscillations decreases 
with decreasing pulse length due to the reduction of the 
effective number of returning electron trajectories. 



IV. INTENSITY DEPENDENCE 

Previous work [l6|, [l?], EH studied the intensity depen- 
dence of the HHG yield at a fixed value of fundamental 
wavelength A c . On the other hand, we have so far fo- 
cussed on the wavelength dependence at a fixed value of 
intensity (1.6 x 10 14 W/cm 2 ). We extend now this anal- 
ysis to the two-dimensional parameter plane (A c ,io) m 
order to explore the underlying mechanisms in more de- 
tail. An example for hydrogen (Fig. for a narrow in- 
terval of wavelength (1 /im < A c < 1.1/mi) and intensity 
(1.3 x 10 14 W/cm 2 <I < 1.6 x 10 14 W/cm 2 ) displays reg- 
ularly shaped ridges each of which can be mapped onto 
a fixed channel closing number R. This regularity is also 
reflected in the cuts through this two-dimensional data 
for different fixed intensities for both hydrogen ( Fig. [5] 
(a)) and argon (Fig. [5] (b)). 

Not only the peak positions but also detailed struc- 
tures of the dependence on A c are quite robust against 
the variation of Jo, when expressed in terms of R. For 
later reference we stress that these remarkable observa- 
tions hold true only when the channel-closing (CC) num- 
ber R is determined with the true ionization potential 
(Eq. (TT])); the use of any other value of effective ioniza- 
tion potential would shift each peak and consequently 
each curve in Fig. [5] by a different amount. This can also 
be understood from the fact that lines of constant values 
of (U p + I p )/huj (with e.g. I p = 10.5 eV) in Fig. (b) 
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FIG. 2: (Color online) Variations of the integrated harmonic 
yield (20 to 50 eV) in a narrow range of A = 1000 — 1100 nm, 
as a function of R, for H. a) comparison between the TDSE 
solutions with the Peaceman-Rachford (PR) and the pseudo- 
spectral (PS) methods, b) the results of the Gaussian model 
with a — 1 a.u. (solid) and 2 a.u. (dotted), c) build-up of 
the interference pattern with increasing number of quantum 
trajectories within the SPA. In b) and c) the vertical axis is 
in arbitrary units. 
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FIG. 3: (Color online) Variations of the integrated harmonic 
yield (20 to 50 eV) in a narrow range of A = 1930 — 1955 
nm, as a function of R, for Ar. a) the TDSE solution, b) 
build-up of the interference pattern with increasing number 
of quantum trajectories within the SPA. 
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FIG. 4: (Color online) Fluctuations of the harmonic yield AY 
as a function of the fundamental wavelength A for hydrogen. 
Solid: 8-cycle flat-top with a 1/2 (1/2) cycle ramp on (off), 
dotted: 14-cycle sin 2 pulse with a FWHM of r p = 7 cycles, 
dash-dotted: 16-cycle Gaussian with a FWHM of r p = 7 cy- 
cles. Other pulse parameter are the same as in Fig. [T] 



deviate from the ridges which manifest as peaks along 
lines of constant Io (Fig. [6] (a)). Results for argon (Fig. [6] 
(b)) show a similar behavior, indicating the applicability 
of the parameter R independent of the atomic species. 



V. MODULATION PERIOD 

The modulation period SX of the harmonic yield is a 
function of the central wavelength A c itself. With in- 
creasing A c , 5X decreases from about 30 nm near 800 nm 
wavelength to w 6 nm near a wavelength of 2 /im (Fig. 
[7]). However, expressed in terms of the channel closing 
number R, the separation of the principal peaks is given 
by 5R = 1 both for the TDSE and the SPA results (see 
Figs. [2] and [3]). Interference peaks appear with this spac- 
ing regardless of intensity (Fig. [5]). The peaks in the 
TDSE results are, however, not located at integer values 
of R, as opposed to the SPA results as well as previous 
theoretical work [It], EBl • This problem was previously 
encountered in the intensity-dependence of HHG and ATI 
[3lT | . In order to recover integer values the use of an effec- 
tive ionization potential I p in place of I p in Eq. ([T]) was 
proposed based on arguments that either the enhance- 
ment was due to multiphoton resonances with pondero- 
motively upshifted Rydberg states [32( (I p corresponds to 
the excitation energy of the resonant state) or that high- 
lying atomic states are strongly distorted by an intense 
laser field to form a quasicontinuum, effectively lower- 
ing the ionization potential 0- As long as one considers 
only the intensity-dependence at a fixed wavelength, the 
difference AI p = I p — I p causes a constant shift of R 
by AIp/oj. Consequently, integer values of R could be 
restored along this axis for a suitable choice of I p . How- 
ever, considering now the wavelength-dependence at a 
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X c [nm] 

FIG. 5: (Color online) TDSE-calculated integrated harmonic 
yield between 20 and 50 eV for H (8-cycle flat-top pulse) in the 
(\ c ,Io) plane. In the contour plot (lower panel), white lines 
show values of constant (U p + I p )/huj, shifted from integer 
values by +0.52, while black lines (only three are shown for 
clarity) represent values of constant (U p + I p )/hw with I p = 
10.5 eV. 

fixed intensity, AI p /uj itself would depend on A. There- 
fore, if the modulation period SX corresponds to 6R = 1 
for the true I p , any other choice of I p different from I p 
cannot shift all the peaks uniformly to integer values of 
R. 

In a further step, we enumerate all the principal peaks 
in Fig.[T](b) fromp = 15 to 117, and plot the CC number 
R p as well as the mismatch to the nearest integer, 

AR P = R P - [R p ], (6) 

of each peak in Fig. [51 The slope of the line fitted to the 
data calculated with true I p (filled circles in Fig. [3a)) 
is nearly equal to unity (w 1.00), while those with I p = 
10.5 eV (diamonds in Fig.[8ja)) have a slope smaller than 
unity. Moreover, although some fluctuation is seen, the 
values of AR p are roughly constant, most of them being 
distributed between 0.3 and 0.6. 

As can be seen from Figs.[2H6l the harmonic yields AY 
are not only composed of peaks separated by SR — 1, but 
also often contain finer structures with sub-peaks. This 
is even more so for longer pulses. In order to extract 
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FIG. 6: (Color online) Wavelength dependence of the inte- 
grated harmonic yield (20 to 50 eV) in the range of A ~ 
1000 — 1100 nm, expressed in terms of R, for a) H and b) Ar, 
for 8-cycle fiat-top pulses for different intensities indicated in 
the figure. 
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FIG. 8: (Color online) Position R p of the principal peaks 
(p = 15 - 117) from Fig. ,: (Up + I P )/huj with the true 

I p = 0.5 a.u. (13.605 eV), A: (U p + f p )/hu with I p = 10.5 
eV. The slope obtained by line fitting is also indicated, (b) 
Corresponding AR P values (Eq. ©) with the true I p = 0.5 



the periodicity of these structures quantitatively we cal- 
culate the power spectrum of AY(R) x A where the 
multiplication by A 5 3 removes the smooth global decay 
(Fig. [1] (b)). We can clearly identify the sharp dominant 
frequency component precisely at Q = 1, corresponding 
to SR = 1. A beat-like structure of a period of ~ 20 seen 
in Fig. Q] (b) gives rise to an additional small side band. 

The Fourier spectrum clearly underscores that the 
peak separation corresponds to SR = 1 throughout the 
entire wavelength range between 800 nm and 2 jj,m. 
This spacing is closely related to the spacing SX in the 
wavelength dependence of the peak positions. Hence, we 
can derive the scaling of SX with A to obtain, 



SX 



const. 



1240 



Jp(eV) + 2.8 x 10- 19 /(W/cni )A 2 (nm) 



nm.(7) 



FIG. 7: (Color online) Variation of the modulation period SX 
with the driver central wavelength A c for atomic hydrogen. 
■ : TDSE, dashed line: SR = 1 (Eq. ©)• 



Equation [7] reproduces the TDSE-calculatcd A depen- 
dence of SX quite well (Fig. [7]). 

The present results as well as those in Sec. II VI strongly 
indicate that the wavelength and intensity dependence 
of the HHG yield calls for an explanation in terms of R 
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FIG. 9: Power spectrum of AY(R) x A ' (see Fig. Q] (b)). Channel-closing number (U p +I p ) / to 



calculated from the true I p in spite of the pronounced 
shift of the peak position from integer R values. 



VI. PEAK SHIFT FROM INTEGER R VALUES 

While the peak separation is given by SR = 1, en- 
hancements do not appear at R = N, with N being 
an integer but shifted by an amount ranging from 0.3 
to 0.6 (see Fig. [8] (b)). This is in clear contrast to the 
SFA prediction for the CC peaks in the literature [l8[ 
as well as to our present SFA results in Figs. [2] (b)(c) 
and [3] (b) . The fundamental difference between the SFA 
and the full solution of the TDSE is - apart from nu- 
merical or analytical solution strategies - that, in the 
former one neglects the excited states and the effect of 
the atomic potential to continuum electrons, which may 
be a serious deficiency for long-range potentials such as 
the Coulomb potential. Recent studies on ionization dy- 
namics and doubly-differential photoelectron momentum 
distributions of hydrogen have shown the significance of 
the long-ranged Coulomb potential in laser-atom interac- 
tion and have illustrated the failure of the SFA near the 
threshold HEUst 

A. Truncated Coulomb potential 

In order to explore the significance of the long tail of 
the Coulomb potential for the present case of interfer- 
ences in the HHG yield, we perform calculations with a 
truncated Coulomb potential, given by 

Veff(r,r c ) = l~ :- ( ,- re)/ r d | r < rc J , (8) 

I r ( r > rc > 

where the effective range of the truncated Coulomb po- 
tential r c is varied between r c = 10 and r c — 70 a.u. 



FIG. 10: (Color online) Comparison of the wavelength- 
dependence 800 nm < A c < 900 nm of the harmonic yield 
AY for H, expressed in terms of R, calculated with the full 
Coulomb potential (marked as "FULL") with the truncated 
Coulomb potentials for varying values of r c as indicated. The 
pulse has a 16-cycle flat-top shape, other pulse parameters 
are the same as in Fig. [T] 

and the width of the cross-over region r<i is chosen to 
be rd = 10 a.u.. For these parameter values, the ion- 
ization potential and the first excitation energy remain 
unchanged to an accuracy of ss 10~ 9 and w 10~ 3 , re- 
spectively. It should be noted that the classical elec- 
tron quiver motion amplitude is a q = 26.3 a.u. for 
I = 1.6 x 10 14 W/cm 2 and A c = 900 nm, and a q = 39.3 
a.u. for A c = 1100 nm. 

We thus explore the entire range from r c /a q « 0.3 to 
r c /ot q ~ 2.3. Convergence to the solution employing the 
full Coulomb potential is reached only for r c as large as 70 
a.u. (see Figs . [TU1 and [TT|) . Most important in the present 
context is a systematic, almost rigid shift of the peaks 
as a function of r c . Only for small r c (sslO a.u.), the 
maxima are found near channel closings (near R equal 
to an integer), in agreement with the SFA results [LSI ]. 
This observation indicates that the Coulomb potential is 
indeed responsible for a - to first approximation - mono- 
tonic and nearly uniform shift of the peaks. It should be 
noted that the long-ranged Coulomb potential manifests 
itself in two seemingly different effects. Firstly, it sup- 
ports high-lying Rydberg states which converge to the 
continuum at threshold. Furthermore, Coulomb scatter- 
ing and deflection influences the motion of the returning 
electron, even at large distances from the core. 

B. Effective ionization potential 

The influence of the potential form on the position 
of the CC was previously identified within the frame- 
work of a 1D-TDSE model [H, It was suggested 
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FIG. 11: (Color online) Same as Fig.QJj] but for a pulse length 
of 8 cycles (flat-top). 

to use an "effective" ionization potential I p in Eq. ([TJ 
when comparing TDSE calculation with models employ- 
ing zero-range potentials to account for high-order above- 
threshold ionization (ATI) spectra at R ^ N [31]. Em- 
ploying an effective ionization potential I p in an SFA 
model roughly leads to a rigid horizontal shift of the in- 
terference structure of the HHG yield, in accordance with 
our observations in Figs. ITD1 and [Til 

Different lines of arguments are invoked for employing 
I p rather than I p . However, they all have in common 
that the existence of a strongly distorted, continuum-like 
excited state e n is considered responsible for an effec- 
tively lower ionization threshold. For convenience, let us 
therefore define AI p = I p — I p , which is expected to be 
a negative quantity (AI p < 0). Different choices of AI p 
are explored. Faria et al. [33j argue that e n should be 
given by the condition that its radius r n « 3n 2 /2 for 
principal quantum number n should match the quiver 
amplitude a q — Fq/u> 2 . Together with the Rydberg en- 
ergy £„ k —Ip/n 2 this would imply 

A/^-^oc/- 1 /^. (9) 

Accordingly, the change of the effective ionization poten- 
tial, AI p , becomes wavelength dependent. On the other 
hand, we have found a fairly rigid equidistancy SR = 1 as 
well as a nearly constant shift AR of the latter away from 
the integers over a wide range of A c . Consequently, if we 
assume that AI p compensates for Ai? p , these quantities 
must satisfy the relation AI p = —(AR p + m)hu>, with m 
being a possible integer offset and AR p defined in Eq. 
^ . Figures [TU] and [TTJ clearly show that the amount of 
the peak shift in R is smaller than unity, hence m = 0. 
This leads to 

AI p — —ARphuj oc A -1 , (10) 



Obviously, hypothesis Eq. ([9]) is not consistent with Eq. 
PU|) . In addition, no upper limit for AI p according to 
Eq. §§§ was discussed in literature. This may lead to the 
obviously incorrect conclusion that AI p ^3.4 eV as soon 
as in a low intensity and low wavelength limit the n = 2 
Rydberg state (or even the ground state!) would govern 
the effective threshold invoked. 

An alternative proposal put forward by Frolov et al. 
relates the energy e n to the formation of an effec- 
tive continuum by broadening of the level with principal 
quantum number n. Accordingly, e n is determined by 
the condition T n = Ae ni where the width T n (related to 
ionization rate) approaches the level spacing Ae„. While 
in the limit of quasi-static tunneling, the tunneling rate 
r„ oc exp I — 2(2|e„|) 3 / 2 /(3F )] strongly depends on the 
field strength F , but only very weakly on A. The result- 
ing value of AIp is estimated to be -3.1 eV for atomic 
hydrogen and I = 1.6 X 10 14 W/cm 2 in Ref. Q- This 
does not meet the condition Eq. (fT0"|) , according to which 
|Aip| should be smaller than the photon energy hw (jl.5 
eV in the present parameter range) of the driving laser 
pulse. 

Moreover, with the help of Eq. (fTD|) we can determined 
the effective parameter dependence of AI p employing the 
numerical values for AR over a wide range of A c and Fq 
(Fig.Q^J). Our results suggest a weak dependence of I p on 
both the wavelength and the intensity, the latter being 
roughly proportional to A -0 3 . This supports neither the 
explanations of Frolov et al. nor of Faria et al. . 

C. Coulomb-corrected classical trajectory model 

In addition to the ability to support (an infinite num- 
ber) excited bound states, the Coulomb potential affects 
the propagation of the rescattering electron which is ne- 
glected in the SFA as well. As the quantum interference 
of electron trajectories is responsible for the oscillation 
in the harmonic yield, their distortion by the potential 
may be crucial. 

In SFA, the time-dependent dipole moment d{t) can 
be expressed as [29( 

<*(*/) = 2 bion ^ ■ e' iSp{u ^^ ■ c ICC (t f ) + ex., (11) 

P{U) 

i.e., a sum over paths P that start at the moment of 
tunnel ionization U with amplitude bi on (ti), evolve in the 
laser field - acquire the phase e - lS p( t i.*/) _ anc l recombinc 
upon rescattering at the core at time t f with the ampli- 
tude c rec (tf). The interference oscillations in the HHG 
yield are controlled by the semiclassical action Sp of the 
path P, which reads: 

S P (U,t f ,I p ) = ^ {jP + A V ))2 dt l + I p (t f -ti), (12) 

where A(t) is the laser vector potential, and p is the clas- 
sical momentum of the returning trajectory. The effect of 
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FIG. 12: (Color online) I p (left axis) and AI P (right axis) as 
obtained from Eq. (|10[) for a broad range of driver wavelength 
and various intensities, a) as function of A c , b) as function of 
intensity near A c = 1950 nm. The solid line compares to the 
prediction of Faria et al. [331 ]. while the dashed line shows a 
power law fit to our data (AI P tx / _0 ' 3 ). 



the Coulomb potential can be incorporated into Eq. (Q2 
with help of an eikonal approximation as a correction to 
the (action) phase with [29[, 



AS P (ti,t f ) 



V EI (r(t'))dt', 



(13) 



Clearly, the eikonal approximation would fail at small 
distances from the nucleus. This difficulty can by by- 
passed using the observation (Sec. lVI A")) that at a cut-off 
r c = 10 a.u. the SFA limit of channel closings at integer 
values of R is reached. Consequently, we set 

V EI {r) = V eff (r,r c = oo) - V eff (r,r c = 10) (14) 

when calculating the long-range phase correction. 

We evaluate Eq. (fT5|) along classical trajectories in the 
laser electric field F(t), confined in the xy-plane, starting 
and ending at the "tunnel exit" zq = I p /Fq. Trajectory 
modifications due to the Coulomb potential are small and 
can be neglected to first approximation (29[, i.e., we use 
the same sets of (ti,tf) as in Eq. (JT2J) . 
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FIG. 13: (Color online) I p (left axis) and AI P (right axis) 
as a function of driver wavelength, employing Eq. [15] and the 
potential of Eq.rj4] Intensity is I = 1.6 x 10 14 W/cm 2 . Lines 
stem from the six shortest orbits recolliding with 20eV, re- 
visiting the core N c (here: zero, one, or two) times before 
recombining (as indicated). For each energy and each N c , a 
short and a long orbit exist. ■: TDSE data (see also Fig.ll2|l. 



It is now suggestive to express this additional phase in 
terms of a change in the "effective ionization potential" , 
I p . Accordingly, 



AS P (ti,t f )/(t f -U). 



(15) 



Figure Q]|] shows I p obtained by Eq. [H] for several (the 
shortest) classical trajectories which contribute to the 
harmonics near 33.6 eV (hence near the center of the 
HHG yield range considered in this work). 

Remarkably, most trajectories (save the shortest one) 
behave qualitatively very similarly. In spite of its over- 
simplification, this Coulomb-corrected model explains 
the behavior of I p even quantitatively well, which is a 
strong indication that the effect of the Coulomb potential 
on the rescattering electronic motion is key to the under- 
standing of the apparent peak shift in the wavelength- 
dependence of the HHG yield. 



VII. CONCLUSIONS 

Using full numerical solutions of the time-dependent 
Schrodinger equation, we have found that the fundamen- 
tal wavelength dependence of HHG with few-cycle pulses 
in the single-atom response features surprisingly strong 
oscillations on fine wavelength scales with modulation pe- 
riods as small as 6 nm in the mid-infrared regime near 
A = 2 /im. Thus, even a slight change in fundamental 
wavelength leads to strong variations in the HHG yield. 
This fine-scale rapid variation is the consequence of the 
interference of several rescattering trajectories with long 
excursion times, confirming the significance of multiple 
returns of the electron wavepacket [lj| . 
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The present oscillations are closely related to similar 
regular peak-like enhancements of harmonic yield as a 
function of intensity Iq [TH, [H, [13, EH, previously dis- 
cussed in connection with channel closings. Our analy- 
sis of the simultaneous wavelength-intensity-dependence 
has revealed that the spacing between adjacent peaks (ex- 
pressed in terms of the channel closing number R) is very 
accurately given by SR = 1 over a wide range of A and 
Iq. This corresponds to the spacing of adjacent channel 
closings as predicted by the strong-field approximation 
(e.g. [ill). The condition SR = 1 holds only if R is de- 
fined with the true ionization potential. However, the 
peak positions are significantly shifted relative to integer 
values. The parametric dependence of the peak shift on 
the wavelength A and the intensity Iq has been investi- 
gated. Our analysis shows that this peak shift can be 
accounted for by the effects of the Coulomb tail on the 
motion of the returning electron. 
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